Electronically Induced Anomaly in LO Phonon Dispersion 
of High - Tc Superconductors 
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The strong, electronically induced anomaly in the spectrum of the longitudinal optical (LO) 
phonons propagating along the main axes of the Cu02 plane is tentatively attributed to the oxygen- 
oxygen charge transfer between the two oxygens in the plane. It is argued that this charge transfer 
can be large and that it is strongly coupled to the zone boundary LO phonons. The corresponding 
negative contribution to the free energy is quartic in the LO phonon amplitude, making the LO 
phonon unstable through the first order phase transition, with a concomitant domain structure. 
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The discovery of the electron-induced anomaly in the 
dispersion of the oxygen LO phonons, which propagate 
along the main axes of the Cu02 plane in La2-2:Sr^Cu04 
[0, P and YBa2Cu3 06+a; has renewed interest in 
the origin of this anomaly and its possible interplay 
with high- Tc superconductivity. Similar anomalies seem 
to occur in La2CaCu04+a; 0, Pb2Sr2(Ca,Y)Cu308 and 
Lii_a;T2_2;04 The oxygen LO phonon is anomalous 
even in the superconducting perovskite Bai_a;K2:Bi03 
1^. On the other hand, structural instabilities of 
a somewhat different nature have been observed in 
La2-^Ba^Cu04 § and Tl2Ba2CaCu208 @, where the 
oxygen displacements occur perpendicularly to the Cu02 
planes. 

The common feature observed in oxygen LO phonons 
of La2-xSr^Cu04 (LSCO) and YBa2Cu306-Hx (YBCO) 
is the occurrence of a break in the dispersion at 
27r/a(0.25,0), half-way between the origin (0,0) and the 
edge (0.5,0) of the Cu02 planar Brillouin zone (BZ). This 
break was tentatively attributed to the (quasi)static 
longitudinal dimerisation of the oxygens with wave vector 
[0.5,0] shown in Fig. la. The latter should produce the 
new BZ edge at [0.25,0] and, correspondingly, a gap in the 
LO dispersion. However, the (quasi)static dimerisation 
has never been observed. This led to the speculation]^ 
that the dimerisation is highly disordered, with various, 
usually very low estimates of the corresponding correla- 
tion length. The length in question should not however 
be shorter than a few times 4a, if an LO phonon gap with 
wavelength 4a is to be attributed to the 2a dimerisation. 

Although the LO phonon anomaly depends strongly 
on doping, the electronic mechanism that produces it 
has not been understood until now. The conventional 
Peierls-like explanation is unlikely, because it requires the 
appearance of the strong [0.5,0] charge density fluctua- 
tions, coupled to the lattice by a linear electron-phonon 
coupling. Although some phonon induced features seem 
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FIG. 1: (a) The (7r,0) oxygen LO phonon. (b) The homo- 
geneous component of the charge fiuctuations induced by the 
(tt, 0) oxygen LO phonons. 



to exist [ [i(|in the electron spectra, recent ARPES mea- 
surements]^ ^ |l^, |l^ have however yielded Fermi sur- 
faces of LSCO and YBCO which do not possess the nest- 
ing properties favoring [0.5,0] charge fluctuations. On the 
other hand, there is no apparent reason why the bare lin- 
ear electron-phonon coupling would itself be particularly 
strong for the LO dimerisation. 

The aim of this note is to examine a mechanism which 
can explain the LO anomaly in LSCO and YBCO. When 
the linear electron-phonon coupling is ruled out, the 
natural next step is to consider a quadratic electron- 
phonon coupling and the charge fluctuations which 
couple through it to the lattice. Such a quadratic 
electron-phonon coupling was already used to explain the 
LTO/LTT transition in La2-a;Ba:ECu04 [|l§, in which 
the oxygens move perpendicularly to the plane. When 
the LO dimerisation mode is involved quadratically in 
the electron-phonon coupling, it can only couple to the 
homogeneous charge redistribution within the unit cell. 
This charge redistribution can be either symmetric or 
antisymmetric with respect to the symmetry operation 
X ^ y. We will show here that the antisymmetric 
charge redistribution is large, that it is exactly the same 
as the one which leads to the LTO/LTT instability in 
La2_2;Baa;Cu04, and that it explains satisfactorily the 
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LO instability in LSCO and YBCO. 

The preceding symmetry analysis can be specified 
within the tight-binding model, which distributes the 
holes over the Cu 3d^2_y2 orbital and the Ox.y "ipx.y or- 
bitals. This model is widely used in the theory of high- Tc 
superconductors either in the small U or in the large U 
limit, where U is the Coulomb repulsion on the Cu site. 
Only four single particle parameters are usually used in 
addition to U, namely, the Cu and Ox,y site energies Ed 
and £p'^, the Cu-0 overlap t and the O^-Oy overlap t'. 



The variation of the crystal-field potential of the ionic 
lattice, introduced by the dimerisation modes with the 
amplitudes Ujr.o and uo.tt, leads in the first place to the 
variation of the site energies, i.e. to a variation of Apd — 
Sd + {£p+£p)/'2 and App = Sp — e^. The variation of the 
oxygen site energies can be calculated in the point-charge 
approximation, and written in the form [nSl E9[ 
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where Q — (7r,0) or {0,7r). Previous calculations 
estimated and from the crystal field potential at 
the oxygen sites. The values of and obtained in 
this way for LSCO and YBCO are given in Table 1. 



TABLE I: (a) The electron-phonon coupling constants and 
/3p for the (vr, 0) deformation, (b) The corresponding con- 
stants calculated for the LTT deformation of the La2Cu04 
lattice. 
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Singling out the homogeneous component of App, 
found that 



it is 
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6 App couples to the homogeneous charge redistribution 

{Up — 11^)^ shown in Fig. lb. 



between Ox and O 



adding a bilinear term to the Hamiltonian Hq, 



H = Ho + n^p6A;p. 



(3) 



In most of the high- Tc superconductors(YBCO is an 
exception) = in the undeformed lattice and SApp 
lifts the Ep = degeneracy of the Cu02 unit cell, in a 
way similar to the Jahn- Teller (JT) effect. Even in the 
presence of the hybridisation t and t' , the JT effect 
leads to large energetic gains. In the conventional JT 
effect the charge transfer between the two levels, whose 
degeneracy is lifted, is a step function of the splitting 



5 App. In general, the effect of hybridisations t and t' is 
however to make this charge transfer analytic 
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where Xm is the exact charge susceptibility for the Hamil- 
tonian (H), i.e. Xpp is related to the exact correlation 
function {rippUpp) for the Hamiltonian Hq of (^. The JT 
effect is now associated with the expectation that Xpp is 
large. Actually, a singular Xpp in (jj) indicates a nonan- 
alytic nature of the charge transfer n^, i.e. recovery of 
the conventional JT effect, with AF linear in A^l,. 
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The free energy variation, corresponding to (IJ), is 



AF = -\xppS^i' 



It should be noted that, according to Refs.(||) and (|^), 
AF of the equation (||) is quartic in the dimerisation M,r,o- 
However, when Xpp is singular, AF is to be taken linear 



(1) in ^A„„, i.e. quadratic in dimerisation. 



Large Xpp is found every time when the Fermi level falls 
close to the two groups of the quasi degenerate electronic 
states respectively associated with the a and b axes of 
the Cu02 plane, which get split by the finite 5 App. This 
is in particular the case when the Fermi level falls close 
to the van Hove singularities at (0.5,0) and (0,0.5) points 
of the two-dimensional Brillouin zone, associated with 
the single electron propagation in the Cu02 planes. The 
single particle picture can be used with some confidence 
in two limits for the Hamiltonian Hq in (^, namely for 
free fermions [U = 0) and for infinite U {U larger than 
Apd, App, t and t'). In this case, for T w K, the real 
part of the generalized susceptibility, given in terms of 
the vertex functions g(k) and ft.(k), reads as 
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Here /[eL(k)] is the Fermi function and eL(k) is the 
energy of the lowest occupied band The case 

g(k) = /i(k) = 1 describes the density of states at the 
Fermi level, np, averaged appropriately over fc^T, and 
.g(k) = h{k) = 9£L(k)/9App the real part of the intra- 
band pp susceptibility, Xpp- Although it has been shown 
that interband effects might also be important in 
high-Tc materials, for the sake of simplicity we shall re- 
tain here only the intraband contribution to Xpp- As- 
suming that A^ = at the outset and for Ci? x Evh, it 
follows than straightforwardly that 



-Xpp « c • ni? « c • log 



max[kBT, {ep - £vh\] ' 



(7) 



with the constant c typically of order 1/5, as illustrated 
in Fig. 2. eg is the energy scale which characterizes the 
van Hove singularity at SyHi assumed two-dimensional, 
logarithmic, rather than extended one-dimensional, i.e. 
more singular than the logarithm. Xpp of Eq. (^ is 
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singular at low temperatures for ep ^ £vH, when the ex- 
pansion (H) breaks down and Ai^ becomes linear in (SA^, 
i.e. quadratic in the dimerisation displacement u. This 
reasoning can be applied to YBCO by the replacement 
(5A^ ^ A^ + (5A^ in Refs.(|) and (|), assuming that 
the value of A^ in the undeformed lattice is sufhciently 
small. 




£F-£p(eV) 



FIG. 2: The density of states and the intraband pp suscep- 
tibihty as functions of the Fermi energy ef, for the simplest 
case, Apd = 0.66 eV, t = 0.73 eV, and = App = (7 = 0. 



when riiT- of (H) is finite, i.e. the negative contribution 
to the free energy AF can be taken as quartic in u. As 
discussed previously , the activation of the negative 
contribution to Ai^ quartic in u occurs through first order 
phase transition. The compelling feature of such expla- 
nation is that the two phases with u = Q and a finite 
u mix, i.e. the (disordered) domains of the dimerised 
phase are expected to condense in the m = phase. This 
might explain the difficulty in the experimental observa- 
tion of the dimerised domains, provided that the latter 
turns out to be small. The theoretical discussion of the 
domain size requires the insight into the gradient terms 
quadratic and/or perhaps quartic in u, which is however 
beyond the scope of the present paper. 

Here, we only wanted to stress that the homogeneous 
charge transfer couples appreciably to the LO dimeri- 
sation because the electron-phonon coupling constants 
of Table I are large and because the corresponding 
susceptibility (0) is large. The disadvantage that the 
electron-phonon coupling is quadratic rather than linear 
in u is set off by satisfactory consequence that the result- 
ing phase transition is of first order, providing a possible 
explanation of the LO anomaly observed in LSCO and 
YBCO. Further experimental and theoretical investiga- 
tions are however required to explain unambiguosly the 
origin of the LO anomaly in the high- Tc superconductors. 
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